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ABSTRACT 

We refine the physical parameters of the transiting hot Jupiter planet XO-lb and its stellar host 
XO-1 using HST NICMOS observations. XO-lb has a radius i?p=1.21±0.03 Rj, and XO-1 has a 
radius i?^,=0.94±0.02 R©, where the uncertainty in the mass of XO-1 dominates the uncertainty of 
Rp and i?*. There are no significant differences in the XO-1 system properties between these broad- 
band NIR observations and previous determinations based upon ground-based optical observations. 
We measure two transit timings from these observations with 9 s and 15 s precision. As a residual 
to a linear ephemeris model, there is a 2.0 a timing difference between the two HST visits that are 
separated by 3 transit events (11.8 days). These two transit timings and additional timings from 
the literature are sufficient to rule out the presence of an Earth mass planet orbiting in 2:1 mean 
motion resonance coplanar with XO-lb. We identify and correct for poorly understood "gain-like" 
variations present in NICMOS time series data. This correction reduces the effective noise in time 
series photometry by a factor of two, for the case of XO-1. 
Subject headings: planetary systems ~ stars: individual (GSC 02041-01657) 



1. INTRODUCTION 

XO-lb is a t ransiting hot Jupiter pl anet orbiting a 
Solar-type sta r dMcCuUough et al.l l2006'). Subsequently, 
iHolman et al.l (|2006[ ) refined the parameters of the star, 
XO-1, and plan et, XO-lb, with ground based opt i- 
cal photometry (iTorres et al.l l2008t1Southworthl 120081 ). 
From a theoretical standpoint XO-lb has the proper- 
ties of a "normal" extrasolar planet, in the sense that 
the measured mass, Mp, and radius, Rp, of XO-lb 
agree with the theoretical expectations after taking into 
account the stellar insolation, and XO-lb likely does 
not contain a subst antial amounts of heavy elements 
(jBurrows et al.ll2007[ ). Spitzer Space Telescope observa- 
tions with IRAC of XO-lb at secondary eclipse detected 
the presence of a temperature inversion in the outer lay- 
ers (1 mba r) of the planetary tem perature-pressure (T- 
P) profile ()Machalek et al.l 120081 ). The wavelength de- 
pendence of the IRAC measurements provides a valu- 
able constraint on the day-side T-P profile and molecular 
abundances. 

Measuring the absorption depth when the planet tran- 
sits its stellar host as a function of wavelength, here 
referred to as transmission spectroscopy, provides ad- 
ditional constraints on the planet's atmosphere. The 
change in transit depth results from opacity variations 
as a function of wavelength in the planet's atmosphere 
and probes the T-P profile in the outermost (< 1 mbar) 
layers of the plane t's atmosphere at the terminator 
(|Fortnev et al.ll2010[ ). To extend transit observations of 

lcburk e@cfa. harvard. edul 

^ Harvard-Smithsonian Center for Astrophysics, 60 Garden 
St Cambridge, MA 02138 

Space Telescope Science Institute, 3700 San Martin Dr., Bal- 
timore, MD 21218 

^ Dept. of Physics and Astronomy, Rice University, 6100 Main 
Street, MS-108, Houston, TX 77005 

* Boston University, Astronomy Dept., 725 Commonwealth 
Ave., Boston, MA 02215 



XO-lb to the NIR, we present high cadence HST NIC- 
MOS observations of XO-lb. The observations were ob- 
tained with the G141 grism, which covers the strong ab- 
sorption features due to H2O, CO2, and CH4. In this 
work, we present new constraints on the XO-1 system by 
summing the fiux over the entire grism spectrum. The 
NICMOS broad-band time series provide a precise mea- 
surement of Rp, stellar density, p^, and transit timings 
for whic h to compare t o prev ious determinations in the 
optical. iTinetti et al.l (|2010f ) present a complementary 
study to this one using the same NICMOS observations 
as analyzed here by enhancing the wavelength res o lution 
at the sacrifice of signal to noise. ITinetti et all polol ) 
examine the relative change in transit depth across the 
G141 grism, and they detect the presence of H2O in the 
atmosphere of XO-lb. 

Currently, NICMOS grism observations for 9 transiting 
extrasolar planets have used >120 HST orbits. Previous 
analyses of grism broad-band time series with NICMOS 
have been limited to a factor of 2-3 t imes worse than the 
Poisso i i expectation (HP 209458b - iGilliland k ArribasI 
(Hods'), 'Swain et al."('2009'); GJ 436b- IPont et al.l (120090 : 
HD 149026b Carter ct al. (2009)). In this work, we iden- 
tify and correct for systematics that improve the relative 
fiux precision for broad-band time series to the expected 
level. The technical details of this work will augment the 
substantial archival NICMOS grism observations and fu- 
ture NICMOS grism observations of transiting extrasolar 
planets. 

The NICMOS observations are described in § H The 
technical details of the data analysis, improvements in 
the treatment of systematics for NICMOS grism obser- 
vations, and a discussion of the expected photometric 
noise level are given in §[3l The XO-1 system properties 
are discussed in § lU and we discuss the results in § [5l 
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The NICMOS grism observations were designed to 
gather high cadence time series of the bright, V=11.2, 
star XO-1 during transits of the planet XO-lb (HST pro- 
gram 10998). Gaps in the time series due to Earth oc- 
cultation necessitate piecing together observations from 
more than one transit. HST observed a transit event on 
Feb. 10, 2008 UT (2454506.6 JD) and again three tran- 
sits later on Feb. 21, 2008 UT (2454518.4 JD). Each visit 
consists of five HST orbits. 

The first-order spectrum of XO-1 from the G141 grism 
(1.1 < A < 1.9 fim) is positioned on the lower left quad- 
rant (quad 1) of the NIC3 NICMOS detector, and it does 
not cross any amplifier boundaries (see Figure [T]) . The 
zeroth-order spectrum lies on the lower right quadrant 
(quad 2), and it enables accurate tracking of the spec- 
trum's position. No other stars contribute significant 
flux in the NIC3 field of view. The Pupil Alignment 
Mechanism was set at —0.53 mm to defocus the Point 
Spread Function (PSF). Defocusing spreads the light over 
more pixels which improves operational efficiency by de- 
laying saturation and imp roves precision by ave raging 
nonuniform pixel response (|Xu fc Mobasheii 120031 ). The 
detector is read out in MULTIACCUM mode with the 
STEP8 sequence and NSAMP=11, resulting in an over- 
aU exposure time 39.953 s for a single MULTIACCUM 
exposure. Including overheads, the exposure cadence is 
49 s between MULTIACCUM exposures. There are ~57 
exposures per HST orbit. A single direct image in the 
F166N filter at the beginning of each HST visit provides 
the reference position of the target for determining the 
wavelength calibration. 

We reduced the data using both custom procedures 
and publicly- available procedures written in IDL. The 
reduction begins from the raw science file (_raw) rather 
than the calibrated science file (_cal). In summary, the 
procedure for reducing an image consists of starting from 
the last-read minus zeroth-read image of a MULTIAC- 
CUM image, applying a wavelength dependent flat fleld, 
and then correcting bad pixels through bicubic spline 
interpolation. These data reduction steps are described 
next, and the resulting images are used to perform broad- 
band photometry in § |3l 

To determine the wavelength dependent flat field and 
to correct for systematic trends in the photometric time 
series (see § [S]) , we map the position of the spectrum on 
the detector. To measure the spectrum location, we fit a 
Gaussian to calculate the centroid of the spatial profile 
for 5 pixel intervals along the dispersion axis. A linear 
fit to the centroids versus position along the dispersion 
axis yields the slope of the spectrum on the detector. 
The lower left panel in Figures [2] & [3] show the spectrum 
slope as a function of image number within an orbit for 
the first and second visit, respectively. Each color rep- 
resents an orbit, and the color code to identify an orbit 
is given in the lower right panel of the figure. The spec- 
trum slope variation within an orbit is much smaller than 
the change between orbits. We measure relative shifts in 
X and y pixel coordinates of the spectrum by cross cor- 
relation of the zeroth-order spectrum with respect to a 
reference image (11-th image in the second orbit of the 
first visit). The zeroth-order position shifts are projected 
to the first-order spectrum's dispersion relation position. 
The projection, Sx = — A cos a and Sy = — A sin a, takes 
into account the spectrum slope, a, and the pixel dis- 



tance between the zeroth-order light and dispersion rela- 
tion's reference position, A = 167.93 pix. The top panels 
in Figure [2] and Figure |3] show the relative x and y pixel 
coordinate shifts of the dispersion relation's reference po- 
sition as a function of the image number within an orbit. 
The dispersion reference position is repositioned within 
0.5 pixel between visits and typically 0.2 pixel within a 
visit. 

We follow t he procedure outlined in 
'Gillil and fc ArribasI (|2003[) to determine the wave- 
length dependent flat fleld, and the wavelength of each 
pi xel is determined usi ng the dispersion relation given 
in iPirzkal et all ()2009D . The centroid of the direct 
image taken at the beginning of each visit determines 
the reference position of the dispersion relation. A flat 
fleld is calculated for each image taking into account 
the reference position shifts and the spectrum slope 
variations. The individual flat flelds are averaged over 
an orbit, and the resulting orbit averaged flat fleld is 
normalized on quad 1 where the first-order spectrum is 
located. The orbit averaged flat fleld is applied to all 
images within that orbit. 

We use four methods to identify the dark, warm, and 
cosmic ray impacted pixels. We begin with the stan - 
dard bad pixel mask for the NIC3 camera (|Sosev II2002I ). 
Second, we identify pixels that significantly vary beyond 
their empirically determined sample standard deviation 
throughout an orbit. Third, warm pixels (> 100 DN) 
in the region for determining the background level are 
identified for correction. Fourth, pixels that deviate by 
> 3(7 from the Poisson and read noise expectation in both 
forward and backward differences are identified for cor- 
rection. On average, each image has 177 pixels (1.1%) 
for correction on quad 1. We replace affected pixels by 
the value of a bicubic spline fitted to surrounding pixels. 

In addition to the spectrum position and slope, we 
make use of the G141 filter wheel telemetry to correct 
for systematic trends in the photometric time series. The 
filter wheel telemetry indicates that the G141 grism posi- 
tion does not return to the same state in between orbits. 
The filter wheel telemetry for each visit (Chris Long, pri- 
vate communication) is given as a function of orbit in the 
lower right panel of Figures [2] & [31 There are two pre- 
ferred telemetry states of the G141 grism (^^0.6 and ~0.3 
fractional telemetry position), and within each of these 
telemetry states the telemetry do not exactly repeat but 
have small variations. In the first visit (Figure [5]), the 
filter wheel telemetry anti-correlates with the spectrum 
slope; The filter wheel telemetry position of 0.6 (orbits 
1 & 3) correspond to low spectrum slope, and the fllter 
wheel telemetry position of 0.3 (orbits 2, 4, & 5) cor- 
responds to the higher spectrum slope. In the second 
visit, the correspondence of fllter wheel telemetry and 
spectrum slope is not as clear. 

3. BROAD-BAND LIGHT CURVE ANALYSIS 

This section describes the process to extract the broad- 
band photometric signal from the calibrated images by 
summing the source flux over the entire flrst-order spec- 
trum. We describe newly identifled systematics of the 
NIC3 camera and the procedure we follow to simultane- 
ously correct the systematics and flt for the transiting 
planet model. The result is an improved estimate for 
the physical properties of the planet XO-lb and star the 
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Fig. 1. — Lower two quadrants of the NIC3 camera showing 
the photometric and background apertures for the target solid and 
blank dashed region, left and right respectively. The image is based 
upon averaging all images in a single orbit. In addition to the 
first-order spectrum in the upper left, the zeroth-order spectrum is 
visible in the upper right. The readout amplifier glow is apparent 
in the lower corners since we analyze raw data rather than the 
CALNICA pipeline calibrated images, which have the amplifier 
glow removed. 
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Fig. 2. — Top Panels: Evolution of the dispersion relation's ref- 
erence X and y position throughout the five orbits of the first visit 
as a function of image number within an orbit. Each color repre- 
sents an individual orbit as given in the key in the bottom right 
panel. Bottom Left: Spectrum slope as a function of observation 
number within an orbit. Bottom Right: Filter wheel telemetry for 
each orbit and color key for each orbit within the first visit. 
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Fig. 3. — Same as Figure [21 but for the second visit. 

XO-1. 

3.1. Photometry and Expected Noise 

We define the source flux to be the sum of counts in 
a rectangular region enclosing the first-order light minus 
an estimate of the background flux contribution. The 
positional shifts of the first-order spectrum (< 1 pix) 



are small enough to adopt the same photometric aper- 
ture and background region for all visits and orbits. The 
photometric aperture is a 113x23 pixel region, and the 
background estimate aperture is 100x65 pixel region lo- 
cated 25 pixels below the photometric aperture. The 
solid rectangular outlines in Figure [1] illustrate the pho- 
tometric and background regions, top left and bottom 
left, respectively. These regions were chosen to minimize 
flux variations in the extracted photometry after correct- 
ing for systematic effects. 

The first-order summed flux from XO-1 is shown in 
Figure IH where the flux relative to the average out of 
transit flux is shown phased at the orbital period of XO- 
Ib. Each HST visit is normalized separately. There is a 
0.6% relative flux difference between the two HST visits 
and smaller, but significant 0.2% relative fiux differences 
between orbits within the same visit. The upper left 
panel of Figure [7] shows in detail the in-transit orbit clos- 
est to orbital phase, (j) = 0.0. For those data, the empiri- 
cal rms noise is tr = 555 ppm per exposure. The expected 
noise from Poisson noise alone is (Tpoi = 180 ppm per ex- 
posure. An additional source of noise arises due to the 
systematic uncertainty in determining the background 
level and read noise. To empirically determine the con- 
tribution of read noise and background level subtraction 
to the photometric precision, we perform the aperture 
photometry on a blank region of quad 2. The target and 
background regions are the same size as for the first-order 
spectrum photometry as illustrated by the dashed line re- 
gions in Figure [1] The expected noise in the blank region 
photometry taking into account the Poisson noise, read 
noise, and the number of pixels in the background and 
photometry regions is 2300 e~ rms. We measure 5400 
e~ rms per exposure variation in the blank region pho- 
tometry, indicating the presence of non-Gaussian noise in 
the background. The empirically determined read noise 
and systematic uncertainty in the determining the back- 
ground level contributes crback=170 ppm noise relative to 
the 3.1x10^ e~ counts from XO-1 per exposure. Adding 
Cback in quadrature with (Tpoi results in the expected rel- 
ative noise for photometry of fitot = 250 ppm per expo- 
sure. The expected noise is 2.2 times smaller than the 
empirically determined noise within an orbit. In addition 
to the higher than expected noise, the residuals qualita- 
tively are distributed uniformly between limits and are 
non-Gaussian in appearance. This broadened distribu- 
tio n of residuals has b een noted previously (see Figure 3 
of (Garter et al.l[20(i9h . 

3.2. Correcting for Gain-like Variations 

The non-Gaussian residuals in the photometric time se- 
ries correlate with 7 preferred states imparted by the de- 
tector electronics that result in gain-like variations. The 
7 states originate from the detector's switching low volt- 
age power supply, which has its clock set to a frequency 
7 times lower than the master clock of the timing pat- 
tern generator (Bergeron, in preparation). The 7 states 
were first identified as additional noise in the detector 
temp erature estimates from the bias level (jPirzkal et alJ 
l2009f ). In this study, we show that the 7 states also im- 
pact photometric time series with NIGMOS. Which of 
the 7 states an exposure belongs to is measurable from 
the zeroth-read image of a MULTIACCUM image. The 
zeroth-read image is obtained immediately following a re- 
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Fig. 4. — Raw relative observed counts for XO-1 phased at the 
XO-lb orbital period. Data is shown from the first HST visit {green 

6 red points) and the second HST visit (yellow & blue points). 
Within a visit, orbits that share the same color share the same 
filter wheel position. 

set of the detector array and it represents the bias level of 
the MULTIACCUM image. For a more accurate measure 
of the relative bias level, the zeroth-read image from the 
first image within an orbit is subtracted, removing any 
bias structure and faint signal of the spectrum recorded 
since the array reset. The mean zeroth-read bias level of 
quad 1 versus the mean zeroth-read bias level of quad 2 
for images within an orbit is shown in the upper right 
panel of Figure [T] The images in this plane form 7 dis- 
tinct groups. The lower left panel of Figure [7] shows the 
same relative photometry as in the panel above it, but 
the measurements that share the same readout state are 
connected by a line. The color coding is the same as in 
the upper right panel of Figure [7l where the 7 states are 
identified. Measurements in the same state systemati- 
cally over or underestimate the relative flux. We correct 
for the 7 state effect by removing the average relative 
flux offset of the state from the relative flux offset av- 
eraged over all measurements in an orbit. The result of 

7 state correction is shown in the lower right panel of 
Figure [71 The rms of the residuals is reduced from 550 
ppm (before) to 288 ppm (after), a factor of 1.9 improve- 
ment. To properly account for the 7 states during times 
when the flux level is varying rapidly (i.e., during ingress 
and egress), the 7 state correction is determined after 
removing the transit model (see below) . 

3.3. Model Fit to Data 

The remaining flux time series, Fotsi are fit to the fol- 
lowing model, 

Fmod 

= i^o X * X $(P,To, Af*,iT!*,5,r/_/y,ii,t2,Wl,U2), 

where i^o x ^ is the correction for the observed flux vari- 
ability due to instrumental effects and intrinsic stellar 
variability between visits (see below) and $ represents 
the relative flux variation due to the transiting planet. 
The orbital period, P = 3.941502 d (Gary, B., pri- 
vate communication), and ephe meris tra nsit midpoint. 
To = 2453887.74679 (HJD; Hol man et al.^ 2 0061. are held 
fixed, but each visit is allowed to have its own transit 
midpoint offset (ti & ^2) from the fixed ephemeris. This 
initial ephemeris is based upon a compilation of profes- 
sional and amateur observations of XO-1. The longer 
time baseline provided by including amateur observa- 



tions allowed a more precise estimate of the ephemeris 
than previously published ephemerides. We adopt a 
qu a.dratic limb dark e ned (u ^ fc U2) transit model as given 
by IMandel fc Agoll (l200l l. and use a G aussian prior 
on th e stellar mass, =1.027±0.06 Mq ([Torres et all 
[20f)8h . The remaining parameters in <& are stellar radius, 
i?^, planet radius to stellar radius ratio, 5 — Rp/ Ri,, and 
first-to- fourth contact transit duration, T/_/y. From the 
constraints pl aced by radial velocity me asurements of the 
XO-lb orbit (IMcCullough et al.l [20061 ) . we assume zero 
eccentricity. The parameter estimates and their uncer- 
tainty are determined f o llowin g the MCMC procedure 
outhned in [Burke et all ([20071 ). The prior adopted for 
Ri,, 5, and tj-iv resu lts in a uniform p rior on i?*, i?p, 
and orbital inclination ([Burke et al.[[2007f) . The priors for 
the other parameters are uniform with upper and lower 
limits well beyond values constrained by the data. The 
exception to this is limb darkening param eters, which 
are ph ysically constrained as described in [Burke et alj 
([20071 ). In particular, we require that the highest sur- 
face brightness be located at the center of the disk (i.e., 
ui > 0). This is consistent with theoretical H-band limb 
da rkening para meters (ui = 0.016 & U2 = 0.441) for XO- 
1 ([Claretll2000f ). The time for each observation is taken 
as the midpoint of the exposure, and is determined from 
the average of the EXPSTART and EXPEND header 
keywords. The resulting modified Julian date is con- 
verted to barycentric Julian date (BJD) on the Terres- 
trial Time system (TT ) using the barycen IDL routing 
([Eastman et al.ir2010f ). The calculated BJD times are 
~1.3 s earlier than heliocentric Julian date times. 

The first orbit of a visit exhibits a systematic ramp up 
in counts over most of the orbit, whi ch likely results from 
telescope or instrument settling ([Gilliland fc Arribas[ 
|2003) ■ Subsequent orbits within a visit also have a ramp 
up effect, but the timescale is much shorter with only the 
first few measurements impacted. To avoid modeling the 
ramp- up effect, we discard the entire first orbit of a visit 
and the first four exposures of each orbit. 

To decorrelate against external parameters to correct 
for the remaining correlated measurements, through trial 
and error, we find that the remaining residuals correlate 
with the external parameters: the spectrum's dispersion 
relation positions, a; & y, and the spectrum slope a (Pig- 
ures[llfc[21). Thus, 

Fo X = Fo X (1 + Co + cix + C2X^ + csu + C4Q;), (2) 

where co,...,C4 are the free linear decorrelation coeffi- 
cients and Fo is arbitrary and fixed (to avoid the degen- 
eracy with Co) at the average out-of-transit flux of both 
visits. We empirically find that a successful correction 
for correlated measurements using the terms in Equa- 
tion [21 requires treating data from the two filter wheel 
states separately (see § [21 for a discussion of the filter 
wheel states). In Figure HI the green and red points rep- 
resent data obtained in the first HST visit, and the yellow 
and blue points represent data obtained in the the second 
HST visit. In each visit, the orbits are assigned to the two 
preferred states of the filter wheel positioning (see lower 
right panel of Figure [2| & [2]) ■ Each visit and filter wheel 

^ http:/ /astro. uni-tucbingen.de /software/idl/aitlib/astro/barvc en. html 
; Input ephemeris file from |http:/ /www. physics. wise. edu/~| 
craigm/idl/down/JPLEPH.405 
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state has its own set of external parameter decorrelation 
coefficients, co,...,C4. Thus, in Figure |4l the orbits that 
share a color share the same decorrelation coefficients, 
however as noted previously, the first orbit for each visit 
is not included in the MCMC analysis. 

There are 5 decorrelation coefficients for each of the 
4 visit/filter wheel combinations, resulting in 20 free pa- 
rameters in 5'. Including the 8 free parameters specifying 
the physical parameters of the system, the overall model 
has a total of 28 free parameters. Before calculating the 
likelihood of a model for a given set of 28 parameters, 
the correction for the 7 readout states is applied. The 
7 readout states are determined for each orbit indepen- 
dently from the residuals of the model. Fobs — Fmod- The 
7 state correction removes the average flux residual for 
each state with respect to the average residual across all 
measurements within an orbit. Thus, the 7 state correc- 
tion does not change the overall flux level of an orbit. The 
likelihood is modeled as independent Gaussian residuals 
with uncertainty, a = 250 ppm. 

The parameters and their uncertainties are based upon 
an overall MCMC chain of 10^ steps after a burn-in pe- 
riod. The scale of the proposal steps in each param- 
eter are set using an automated iterative algorithm of 
proposal step size adjustments until the acceptance frac- 
tion, 0.2 < / < 0.3, is reached. After an initial burn-in 
period to determine the proposal steps, the proposal dis- 
tributions are flnalized, and the results are based upon 
the remaining steps. The longest auto correlation length 
amongst the parameters is 600 steps. Initial t ests with 
a para llel tempering (PT) MCMC algorithm (|Gregorvl 
l2005al) with 7 parallel chains, did not show any evidence 
for multimodality amongst the 28 free parameters. Thus, 
the more time-consuming PT MCMC was not needed to 
explore the parameter space and helps verify that the 
single chain reliably explored the parameter space. 

4. RESULTS 

The resulting transit light curve, residuals, and cor- 
rections are shown in Figure [51 The top panel shows 
Fobs/'^, the relative flux of the observations after divid- 
ing out the best fit (in a sense) decorrelation function 
with correction for the 7 states. The color coding of the 
points is the same as in Figure 21 The middle panel 
of Figure [H shows the residual relative flux of the ob- 
servations from the complete model. The resulting rms 
residual, Crms = 240 ppm, is slightly less than the ex- 
pected uncertainty, Ctot = 250 ppm (see § [3]) . We pre- 
sume the latter results from a slight overestimate in the 
contribution of cTback to crtot- The lower panel of Figure [H 
shows the relative flux correction due to The peak to 
trough relative flux variation in ^ reaches 0.16% for or- 
bits within the same visit and fllter wheel position. For 
display purposes, ^' is normalized to its average value 
for each visit/filter wheel combination. Within a visit, 
the external parameter decorrelation normalization, cq, 
varies by 0.12% in relative flux. In the 12 days between 
the two visits, the observed flux of XO-1 varied by 0.6% 
due to intrinsic variability of the star or the instrument, 
but we cannot distinguish between these two possibili- 
ties from these data. The minimum = 362.8 with 362 
degrees of freedom indicates the model is an acceptable 
flt to the data. 

Tables [1] & [2l show the resulting parameter estimates 



and their uncertainty for XO-lb and XO-1 along with 
previous determinations from the literature, respectively. 
The parameter estimates come from the median of 
MCMC samples, and the uncertainties inscribe 68.3% 
of the MCMC samples. The uncertainties include the 
impact of the assumed prior on M^, =1.027±0.06. This 
affects the physical properties of the system that are not 
directly constrained by the light curve (e.g., i?* & Rp)- 

4.1. Light Curve Quality 

The light curve quality is high enough that the uncer- 
tainty in Ri, and Rp is dominated by the uncertainty in 
Mi^. A solution assuming flxed M^, results in cr^, = 0.008 
Rq and tTp = 0.013 Rj uncertainty in Ri, and Rp, re- 
spectively. Table [31 illustrates the sensitivity of the pa- 
rameters that are directly constrained by the light curve 
(i.e., independent of Mj,) to fixing the limb darkening pa- 
rameters at their theoretical expectation. We ad opt the 
H-ban d values (ui = 0.016 & ua = 0.441) from iClaretl 
(2000*). Fixing the limb darkening parameters results in 
a Ax^ = 1.3 worse fit, a non-significant difference in 
the quality of fit. However, the estimate for Rp/R^ is 
la smaller and the uncertainty is 40% smaller. Fixing 
the limb darkening parameters can lead to more precise 
model fits, but in the case of high quality data, it may 
result in lower accuracy when the stellar bri ghtness pro- 
file differs from th e theoretical expectation (jSouthworthl 
l2008HClaretj|2009l ). Also shown in Table[3]is the precision 
possible if the observations were free from systematics. 
The results were obtained by fixing the decorrelation co- 
efficients, Co, ...,04, at the values that minimize x^- The 
need to correct for systematics reduces the precision in 
Rp/ Ri, by a factor of 5 and reduces the precision in the 
transit midpoint by a factor of 1.2. 

The results shown in Table|31emphasize the difficulty in 
judging the light curve quality based upon comparing the 
resulting uncertainty in the model parameter estimates 
alone. Assuming a more precise estimate of Mi, (or equiv- 
alently i?*), fixing the limb darkening parameters, or as- 
suming fixed decorrelation coefficients for in-transit data 
based upon their solution from out-of-transit data, will 
all result in higher precision for the system parameters, 
but such procedures may result in less accuracy. Com- 
paring the quality of light curves should not be done out 
of context. For example, with comparable quality data, 
the precision for a transit model analysis using fixed limb 
darkening coefficients should not be directly compared to 
the precision that results from a transit modeling analy- 
sis that allows the limb darkening coefficients to vary. 

4.2. Understanding NICMOS Systematics 

Recently, iTinetti et all ()2010f ) analyzed the same ob- 
servations of XO-1 to measure the transmission spec- 
trum of the atmosphere of XO-lb. Although their 
analysis focused on individual spectral channels across 
the G141 grism, their independent analysis techniques 
share some similarities with the approach adopted in 
this study. The 7 state correction described in this 
study is closely related t o their correction f or "channel- 
to-channel" correlations () Swain "erall[2008l) . The gam- 
like variations associated with the 7 states, which ap- 
pear to be wavelength independent, would coherently af- 
fect the relative flux residuals averaged over all wave- 
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length channels, and be removed through the "channel- 
to-channe l" corrections . The "channel-to-channel" cor- 
rection of iSwain et al.l (|2008D cannot be applied to the 
broad-band photometry, since by definition the broad- 
band photometry h as a single channel. The study by 
iTinetti et aD ()2010[ ) analyzed only the second HST visit 
to XO-1 (yellow and blue points in Figure 0]), as the first 
visit (green and red points in Figure |4]) was deemed too 
photometrically unstable. By treating orbits in different 
filter wheel positions separately, we are able to provide 
reliable model fits (reduced ^ 1) using data from 
both HST visits to constrain the properties of the XO- 
1 system. Analysis of more NICMOS observations are 
warranted, but it indicates our methodology enables a 
coherent procedure to help analyze NICMOS datasets of 
transiting planets that were previously thought to be too 
photometrically unstable to provide useful results. 

4.3. Impact of Stellar Spots 

The light curve analysis assumes the stellar surface is 
described by the limb darkening function. However, the 
presence of dark spots or bright facula e on the surface o f 
the star violate this assumption (e.g., iPont et aLll2008[ ). 
The discovery photometry (0.8% precision) and slow ro- 
tation (vsin i < 3km s~^) in dicate that XO-1 is not an 
active star (|McCullough et a l. 2006). There is a 0.6% 
apparent fiux difference in the measured counts between 
HST visits; XO-1 appears brighter during visit 2. The 
photor netric stability of the NICMOS cameras is of order 
- 1% (jThatte et al.il2009[ ). and orbit to orbit differences 
of 0.2% in flux within a single visit are present. Thus, it 
is not clear whether the flux difference between the first 
and second visit is due to intrinsic variability in XO-1 
or due to instrumental photometric instability. If the 
r/ = 0.006 decrease in fiux results from the appearance of 
an unocculted dark spot, then the Rp/Ri, as measured 
for the transit when the dark spot is present needs to be 
reduced by a factor of y/1 — r] — 0.997 to compare to the 
measured Rp/R^, when the spot was absent (i.e., the tran- 
sit depth is deeper when more unocculted dark spots are 
present). The expected 0.25% change in i?p/i?* between 
visits due to the presence of a hypothetical dark spot, 
is smaller than our precision with which we currently 
measure Rp/R^,. To verify this expectation, we added a 
free parameter to our model allowing each visit to have 
its own Rp/Ri,. The resulting model negligibly improves 
the fit, Ax^ < 1, however ARp/R^ = -0.0003 ± 0.0009, 
or 0.23% lower Rp/Ri, for the second visit. This agrees 
in size and direction as the expectation, however it is 
within the statistical uncertainty, and thus we don't for- 
mally adjust the results in Table [T] 

4.4. Impact of Correlated Data 

The likelihood used in the Bayesian posterior assumes 
the residuals are independent and can result in underesti- 
mated para meter uncertainties if the residuals are n ot in- 
dependent ()Pont et al.l 12001 ICarter fc WinnI l2009l) . We 
analyze the sample autocorrelation function ( Box et aP 
|2008| ) to verify that the correction for systematics re- 
moves the temporal correlations present in the raw data, 
and that our use of the likelihood that assumes inde- 
pendence is valid (Figure [5]). To reliably determine the 
sample autocorrelation especially for large lags, the sam- 
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Fig. 5. — Sample autocorrelation of the model fit residuals as 
a function of image sequence lag. Upper and lower smooth curves 
indicate the 3cr limits for the null hypothesis that the residuals are 
independent for the sample autocorrelation statistic. 

pie autocorrelation is calculated on the vector of residu- 
als that are ordered by image number rather than time. 
The smooth upper and lower curves in Figure[S]shows the 
3(T limits as to the expectation of the sample autocorre- 
lation statistic if the underlying population of residuals 
are independent (Kan & Wang 2010). The sample auto- 
correlation quickly dies off by lag 1, ai = —0.067, but a 
few of the lags spike above the 3a expectation. 

We further explore the assumption of independence in 
the likelihood through Monte Carlo simulation of the raw 
data and subsequent retrieval of the system parameters. 
We generate a red noise vector that has the same vari- 
ance and autocorrelation function as the model residuals. 
To ensure the sample autocorrelation function does not 
underestimate the true autocorrelation function, the sim- 
ulated red noise has an autocorrelation for all lags (> 1) 
1.5 times the sample autocorrelation shown in Figure [5j 
In addition to the transit model with red noise, we add 
correlations with external parameters and 7 state offsets. 
By comparing the scatter of parameter estimates derived 
from a minimization of the simulated light curves and 
compare to the uncertainty estimates of an MCMC anal- 
ysis of the actual data for fixed stellar mass, we find the 
uncertainties agree within 5% for all the transit model 
parameters except for the uncertainty in R*, which is 
18% larger. However, the larger uncertainty in is 
negligible compared to the much larger uncertainty in 
R^ that results from the uncertainty in M^,. These tests 
show that the use of a likelihood assuming independence 
does not result in a significant underestimate in the pa- 
rameter uncertainties. 

4.5. Transit Timing Variations 

The exposure timing information in the NICMOS im- 
age headers is on a UTC system, which includes leap sec- 
onds. To place timings on a system free of leap seconds, 
we provide the transit timing measurements from this 
study along with previous transit timing measurements 
from the literature on the uniform Barycentric Julian 
Date Terrestrial Time (TT) system in Table HI To re- 
fine the ephemeris of XO-lb, we analyze precise transit 
timing measurements from the refereed literature based 
upon a complete transit and not impacted by strong 
tren ds. A subset of the transit timing measurements 
from lHolman et ahl (|2006( ) and iCaceres et al., (,2009. ) meet 
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these characteristics. Specificahy, the first 6 transit tim- 
ing measurements hsted in Table |4] are included. The 
resulting ephemeris is given in Table [1] and the timing 
residuals from the linear ephemeris model are shown in 
the top panel of Figure [6l With the published uncer- 
tainties on the transit timing measurements, the linear 
ephemeris model is a moderately poor fit to the transit 
times. The resulting = 16.08 with f — A degrees of 
freedom has a 2.98 a chance of occurring randomly. A 
second set of uncertainties for the ephemeris of XO-lb is 
given in Table [T] that result from scaling the published 
transit time uncertainties to enforce a reduced — ^ 
for a linear ephemeris. The significance of rejecting the 
linear ephemeris depends most on transit event 86 in Fig- 
ure |6l Without transit event 86, the linear ephemeris fit 
improves to — 6.94 with v = 3, which has a 1.79 a of 
chance occurrence indicating an acceptable fit. 

Additional transit timing measurements of XO-lb, 
with lower precision, are consistent with the lin- 
ear ephemeris determined from the higher precision 
transit timings. We select additional transit timing 
measurem ents of XO-l b from the Exoplanet Transit 
Database (jPoddanv et al. 2010) and the Amateur Ex- 
oplanet Archive (AXAj^ that have observations based 
upon > half of the transit and not impacted by strong 
trends. We list the additional transit timings that sat- 
isfy these characteristics in Table |4l The transit timings 
are converted to the BJD (TT) ti ming system, and w e 
also include the measurement from iRaetz et al.1 (j2009t ). 
The bottom panel in Figure [6] shows the residual timings 
(Black points) from the linear ephemeris model, which 
was determined from the higher quality transit timings 
(Red points). These additional transit timings have a 
= 23.2 for v — 23 indicating the linear ephemeris 
model is adequate for these data. 

The transit timing residuals from a linear ephemeris 
shown in Figure IHl qualitatively give the appearance for 
a sinusoidal pattern (Gary, B., private communication). 
Adding a long period sinusoidal component to a linear 
ephemeris results in a Ax^ = 11.6 improvement from 
the linear ephemeris, with a period for the sinusoidal 
component of 118.3 transit cycles (463 day). Bayesian 
evidence is one method to quantify the most probable 
model that describes the data. The details to calcu- 
late the Bayesian odds ratio between a linear ephemeris 
model and the linear -|-sinusoid ephemeris is given in Ap- 
pendix |Al Using the transit timing measurements in Ta- 
ble m the odds ratio of the linear -|- sinusoidal model to 
linear model is 5:1. This is substantial evidence in favor 
of the linear -|-sinusoidal model, but much less than the 
100:1 odds ratio typically accepted as a decisive result. 

Overall, the linear -|-sinusoidal model cannot be dis- 
tinguished from the linear model, and given the cur- 
rent data both models are acceptable descriptions of 
the data. The most likely value for Pgin = 6.84 transit 
events (26.77 day), with the previously identified longer 
period Pgin = 120 transit events (470 day), the second 
most likely solution. The preference for a linear -|-sinusoid 
model is weakened, odds ratio 2:1, when the transit tim- 
ing measurement from the top panel of Figure |6] that 
most significantly departs from the linear ephemeris is 

® http://brucegary.nct/AXA/x.htm 
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Fig. 6. — Top: Deviations of high precision transit timings of 
XO-lb from the literature and new timings from this work from 
a linear ephemeris model. From left to right, measurements with 
Keplercam on the FLWO 1. 2m ijHo lman et al. 2006), SOFI on the 
NTT and ISAAC on UTl HCaceres et al.ll200a '). and results from 
the two HST visits reported in this work. The ephemeris of XO-lb 
is calculated from these transit timings. Bottom: Deviations of ad- 
ditional transit timing measurements of XO-lb from the ephemeris 
derived in this work (black points). The timing residuals from the 
top panel are reproduced in this panel {red points). 

not included in the analysis. 

One potential source of Transit Timing Variations 
(TTV) is the dyn amical influence of additional plan- 
ets o rbiting XO-1 (iHolman fc Murray! [2005t lAgol et al.l 
120051 ). In particular, the amplitude of TTV for Earth 
mass planets orbiting at the 2:1 interior or exterior mean- 
motion resonanc e with a hot Jupiter planet like XO-lb 
can reach >50 s (jHaghighipour et al.l 20091. A full anal- 
ysis on the limits for companions to XO-lb from the null 
detect ion of TTV for XO-lb is beyond the s cope of this 
work ()Agol et al.l 120051: lAgol fc Steffenll2007[ ). However, 
we use the Bayesian evidence calculation to rule out the 
presence of Earth mass planets in an exterior 2:1 mean 
motion resonance with XO-lb. We model the TTV of 
XO-lb due to a coplanar companion at the exterior 2:1 
mean motion resonance by adding a sinus oid with an ex- 
pected period of 72. 2P (i Agol et al.ll2005D to the transit 
timing measurements in Table |4l Amplitudes >45 s for 
the sinusoid TTV result in an odds ratio >100:1 for a 
decisive detection. 

5. DISCUSSION & CONCLUSION 

Without the interference of the Earth's atmosphere, 
space based observations provide an opportunity to 
achieve high precision, Poisson-limited photometric ob- 
servations. However, the thermal forcing from the or- 
bital cycle of HST and uncorrected instrumental sys- 
tematics present a challenge to Poisson limited perfor- 
mance. We correct for gain-like variations in high ca- 
dence, high precision time series observations in NIC- 
MOS, which improve by a factor of two the relative flux 
measurement. Although the baseline of observations is 
limited, there is some evidence that the non-repeatability 
in the G141 grism positioning between orbits affects the 
system throughput at the 0.12% level in relative flux. 
When observations taken with similar grism positioning 
are grouped together, only a few terms related to the 
PSF positioning and shape are needed to decorrelate the 
systematics that appear on the HST orbital time scale. 
The remaining source of noise is consistent with Poisson 
and the uncertainty in determining the background level. 
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Despite the improved modeling of the systematics, the 
physical transiting planet model contains some degener- 
acy with the corrections for systematics. The largest im- 
pact is for determining the transit depth, Rp/R^,, which 
has an uncertainty 5 times larger than if the observa- 
tions were of equivalent quality but free from red-noise 
systematics. Parameters such as a/i?* and transit tim- 
ing are le ss sensitive t o the treatment of red-noise sys- 
tematics. iSouthworthI (|2008f ) reaches similar conclusions 
for the analysis of the high quality ligh t curves of HP 
209458b taken with t he STIS on HST (Br own et al1l200TI : 
iKnutson et al.l [20071 ). This emphasizes the necessity of 
simultaneous fitting of the transit model and systemat- 
ics in order to provide a more realistic assessment of the 
uncertainties. 

Despite the limitations from needing to correct for the 
systematics, we greatly improve upon the precision for 
determining the properties of XO-lb directly measurable 
from the light curve such as a/i?^, Rp/Ri,, p*, and transit 
timing. We find no significant difference in the determi- 
nation of Rp/Ri, betwee n the NIR and previous deter- 
minations in the optical (jTorres et al.' 2008; ' SouthworthI 
|2008, Reproduced in Tables [I] fc El). Using the same 
data set as this studv. iTinetti et aD (j2010) independently 
find a similar absorption depth averaged over the full 
wavelength coverage. Using transmission spectroscopy, 
precise and accurate measurements of the transit depth 
as a function of wavelength are sensitive to opacity 
variations in the upper levels (1 mbar) of the plane- 
tary atmosphere, and theoretical planetary atmosphere 
models can be used to interpret these measurements to 
constrain the abu ndance of molecules and temperature- 
pressure profiles (iSeager fc Sasselovl 120001: [Bro wn 2001; 
iMadhusudhan fc Seageril2009l : iFortnev et al.ll2010.) . For 
example, in their anal y sis of G141 grism data for HD 
149026b, iCarter et all (|2009 ^ find a larger absorption 
depth averaged over the G141 grism when compared to 
optical and Mid-IR measurements than the theoretical 
models predict. The technical improvements to NIC- 
MOS data analysis that we outline in this study may 
lead to a more precise estimate of the absorption depth 
in the large archive of NICMOS grism observations of 
transiting extrasolar planets. 

The transit light curve alone does not allow measuring 
Rp, and with the quality of light curve from this work, 
the uncertainty in Rp is dominated by the adopt ed esti- 
mate for Mi, based upon stellar isochrones ( Torre i~et al.l 
120081 ). Alternatively, measuring the parallax of XO- 
1, along with Te//, apparent magnitude, and bolomet- 
ric correction determinations, yields a constraint on Ri,. 
This constraint on R^, and the transit light curve will pro- 
vide a separate estimate for M-^. As part of the proposal 
for these NICMOS observations. Fine Guidance Sensor 
(FGS) on HST observations were obtained to measure 
the parallax of XO-1. We estimate the FGS distance will 
provide an estimate of Ri, with an uncertainty of 5%, 
and this will provide a 15% constraint on Mi, using the 
transit light curve from this study. Asteroseismology pro- 
vides an additional estimate of and has recently been 



demonstrated to yield very precise estimates for the stel- 
lar host and planet in the case of HD 17156b (Nutzman, 
in preparation), and will be a routine procedure for the 
transiting planets f ound around the bri ghter stars with 
the Kepler mission (iGilliland et al.l[2010l ). Measuring the 
bulk mass and radius of planets with the highest preci- 
sion will require improvements in understanding the mass 
and radius of their stellar hosts (|Southworthll2Q08l ). 

The improved precision of the NICMOS time series 
photometry has provided two additional measurements 
of the mid-transit time of arrival. The timing between 
the two HST visits analyzed in this work is different by 
2ct from a linea r ephemeris. Also, a single transit mea- 
surement from (jCaceres et al.|[2009l ) departs most signif- 
icantly from a linear ephemeris. Based upon the current 
measurements of XO-lb transit timings, any bona fide 
TTV for XO-lb will likely have a peak-to-trough ampli- 
tude <90 s, and if the TTV are concentrated in single 
events, then single timing precision needs to be <23 s. 
The HST observations presented in this study achieved 
(7=8.6 s and (T=15 s precisions in transit timing. The non 
continuous nature of HST observations negatively im- 
pacts the transit timing. The timing precision with HST 
can be improved by optimizing coverage of the ingress 
and egress portion of the light curve. However, transit 
timing measurements with (t=13 s precision are possible 
with g round based o bservations as demonstrated on XO- 
lb bv iH olman et al . (2006). More detailed analysis is 
warranted, but the current transit timing data rule out 
coplanar Earth mass companions to XO-lb orbiting in 
2:1 mean motion resonance. There have been no addi- 
tional radial velo city measurements publis hed since the 
discovery paper (jMcCullough et all I2OO60 . Additional 
radial velocity measurements can aid in constraining the 
presence of additional planets orbiting XO-1. 

The Bayesian evidence calculation that we outline 
in Appendix |^ can be applied to other transit tim- 
ing data sets. The simplified modeling of TTV from 
a linear ephemeris with a sinusoidal component makes 
the model selection process numerically expedient. As 
the quantity and quality of transit timing measurements 
expands from ground-based observations and especially 
with the contributions of the CoRo T and Kepl er missions 
(jCsizmadia et al.|[2Q10l : IGilliland et al.ll2010bD . early dis- 
covery of TTV will ai d in priorit izing and planning the 
followup observations (|Fordll2008l ). 

This work benefited from discussions with Philip Nutz- 
man, Dan Fabrycky, Joe Hora, David Charbonneau, and 
Jeff Stys. We thank Mark Swain and Pieter Dieroo for 
conversations about spectrophotometry with NICMOS. 
We thank Josh Winn and Valentin Ivanov for discussing 
their published transit timings, and the following ob- 
servers for their generosity in making their transit timing 
measurements publicly available, Anthony Ayiomamitis, 
Cindy Foote, Bruce Gary, Joao Gregorio, Bill Norby, 
Gregor Srdoc, Jaroslav Trnka, and Tonny Vanmunster. 
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APPENDIX 

BAYESIAN EVIDENCE FOR SINUSOIDAL TRANSIT TIMING VARIATIONS 

Bayesian evidence is one method to quantify t he most probable mod el that describes the data, and we follow the 
procedure as outlined in lGregorvl (|2005al lbl. l2007f ): lFord fc Gregory! ()2007f ) to calculate the Bayesian odds ratio between 
the linear ephemeris model and the linear+sinusoid ephemeris model. The algorithm presented by Gregory (2005a) is 
tailored to quantify the odds for a radial velocity data set to be modeled by no planet, one planet, or > one planet, 
and we independently implement the algorithm for the model comparison problem at hand. We model the transit 
timing ephemeris as 

To + ExPnn + asm(^^ + 4^ , (Al) 

where the free parameters are the ephemeris zeropoint. To, linear ephemeris period, Pun, amplitude of the sinusoidal 
component, a, sinusoidal component period, Psin, and sinusoidal component phase, 0. The independent variable in 
the model is the transit event, E. Thus, we calculate the odds ratio between the linear ephemeris model with only two 
free parameters, To & Pun, and the full linear+sinusoid model with five free parameters as shown in Equation lAll We 
employ the independent Gaussian model for the residuals in the likelihood. We include an additional free parameter 
in the likelihood that scales the reported errors on the transit timing measurements, s. Allowing s to vary enables 
the odds ratio between the competing models to take into account our incomplete knowledge as to the overall scaling 
of the transit timing uncertain t ies. H owever, we are assuming that the transit timing uncertainties are correct in a 
relative sense. iFord fc HolmanI ()2007|) employ the same linear+sinusoid ephemeris model to characterize and simulate 
the detection of the TTV signature from a Trojan companion using a lomb-scargle periodogram approach. 

The priors for To and Pun are uniform since they are already well determined, and their prior range accommodates 
the constraints placed by the data. We adopt a Jeffreys prior for s that is constrained between 0.25< s <4.0. We 
adopt a modified Jeffreys prior for a with a break at 1.0 s and maximum of IC* s. The maximum a is th e largest 
expected amplitude of transit t i ming variations expected in the planetary regime ()Holman fc Murravl [2001 . Based 
upon the discussion in iGregorvl (|2007[ ). we choose a Jeffreys prior for Pgin with a minimum of 0.9 transit events and 
a maximum of 1000 yr. The prior for < < 27r is uniform. The PT MCMC algorithm has 7 parallel chains, 
and an automated iterative algorithm of proposal step size adjustments until the acceptance fraction, 0.2 < / < 
0.3, is reached. The calculation is based upon the MCMC steps after the proposal step sizes are finalized. The 
linear+sinusoid ephemeris model is not optimized for the complicated signature TTV signals that can be present 
and will underestimate the importance of a TT V signature that is concentrated in few events. H owever, a variety of 
TTV signatures are qualitatively sinusoidal fe.g.. lFord fc Holmanll2007l : iHaghighipour et al.ll2009[ ). and in general, the 
dynam ical variations of the osculating elements in resonant interactions are well described by expansions of oscillating 
terms (jMurrav fc Derm ott 1999). 

The Bayesian Information Criterion (BIG) is an approximation to the Bayesian evidence that is most accurate in the 
limit of a large number of measurements with a simple likelihood surface (Kass & Raftery 1995). In this example, the 
number of transit timing measurements is small, and the multimodal likelihood surface invalidates the Taylor series 
expansion approximation implicit in BIG. However, the ABIC=5.2, between the linear and linear+ sinusoidal model 
mildl y prefers the linear+sinusoidal model. ABIG>10 is necessary for a decisive decision with the BIG ()Kass fc Raftervl 
|1995() . Both the Bayesian evidence and BIG show that there is a mild preference for the linear+sinusoidal model, but 
not a definitive one. 
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Fig. 7. — Top Left: Raw relative flux for the single orbit fully in transit as a function of image number within an orbit. The noise is 
2.2 times higher than expected and broadly distributed (i.e., non-Gaussian). Top Right: Zeroth-read bias level from the lower left detector 
quadrant (quad 1) versus the the lower right detector quadrant (quad 2) for the same orbit as shown in top left figure. The points are color 
coded into 7 distinct state groupings. Lower Left: The relative flux observations that share the same 7 state grouping are connected by 
lines with the same color coding as the top right figure. The observations for each state are systematically above or below the mean flux 
level. Lower Right: Each group is forced to the average flux level of an orbit resulting in nearly a factor of 2 improvement in the noise. 
The remaining correlated red-noise is removed by decorrelation with external parameters. 



TABLE 1 
Planet Properties - XO-lB 



Parameter 


This Work 


Torres ct al. (2008) 


Southworth (2008, 20()9) 


Mass, Mp [Mj] 
Radius, Rp [Rj] 


0.92±0.08 
1.21±0.03 


0.92±0.08 
1.21±0.04 


0.94±0.07 
1.22±0.07 


Velocity semiamplitude, K [m s~^] 
Semimajor axis, a [A.U.] 


0.049±0.001 


0.0493±0.0009 


0.0499±0.0008 


Orbital inclination, i [deg] 


88.8±0.2 


88.8±+{]-J 


89.1±0.8 


Scaled semimajor axis, a/R^ 
Planet-to-star radius ratio, Rp/Ri, 


11.24±0.09 
0.1320±0.0005 


11 ^^++0.03 

0.1326±0.0004 


0.1317±0.0019 


Transit duration (l''*-4*'' contact), t/_/^ [hr] 


2.971±0.006 






Ingress/Egress duration (l'''-2"'* contact), ti-u [hr] 


0.366±0.007 






Impact parameter, b 
Planet-to-orbit radius ratio, Rp/a 


0.23±0.04 
0.0117±0.00012 


0.01166±0.00035 


Star-to-orbit radius ratio, R-kja 


0.0890±0.0007 




0.0886±0.0019 


Summed radius-to-orbit ratio, (ij* + Rp)/a 


0.1007±0.0008 




0.1003±0.0022 


Planet gravity, loggp [cgs] 


3.19±0.03 


3.21±0.04 


3.199±0.04 


Planet density, pp [g cm~^] 
Safronov, 


0.64±0.05 
0.073±0.006 


0.65±0.09 
0.0744±0.006 


0.69±0.08 



Orbital Period, P [day] 3.94150685 ± 0.00000091 ± 0.0000018=» 

Time of midtransit, To [BJD(TT)] 2453887.74774 ± 0.00011 ± 0.00022=''*' 



The second set of uncertainties correspond to enforcing reduced = 1 for a linear ephemeris model through scaling the 
uncertainties in the transit timing measurements. 

Time is on th e BJD Terrestrial Time (T T) system. For times on the BJD (UTC) system subtract 65.184 s. 
" Adopted from lMcCuUoueh et all 12009 1 
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Fig. 8. — Top: Best fit relative flux of the observations after correction for the 7 readout states and decorrelation with external parameters. 
Middle: Residual from the best fit transit model and corrections. Bottom: Best fit correction for the external parameter decorrelation. 
The color coding is the same as Figure |4] 



TABLE 2 
Stellar Properties 



XO-1 



Parameter 



This Work 



Mass, M* \Mq\ 
Radius, ij* [Rq\ 
Effective temperature, T^ff [K] 
Metallicity [Fe/H] 
Stellar gravity, logg [cgs] 
Stellar density, [g cm~^] 
First limb darkening coeff., u\ 
Second limb darkening coeff., U2 



1.027±0.06'' 
0.94±0.02 

0.02±0.08='''" 
4.50±0.01 
1.73±0.04 

<i<T[/.L.0.07 
0.35±0.08 



Torres et al. (2008) Southworth f2009) 
1.066±0.051 



0.934±0.035 



4.509±in^? 
1.877±1°:°}5 



0.95±0.025 



4.510±0.018 
1.75±0.11 



Adopted from Torres et al.' C200g) 
^ Adopted from .McCuUouah ct all l|2006l ) 



TABLE 3 
Precision of Light Curve 



Parameter 


limb free 


limb fixed 


Sys fixed 




11.24±0.09 


11.28±0.09 


11.26±0.07 




0.1320±0.0005 


0.1316±0.0003 


0.1316±0.0001 


Ti-iv [dy] 


0.1238±0.0003 


0.1240±0.0002 


0.1240±0.0002 


ti [s] 


(T=8.6 


cr=8.6 


(T=6.9 
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TABLE 4 
Transit Timing Measurements 



BJDo (TT)=' 


a 




Source 




0.00010 
0.00017 


_L lllo VVUi K 
± lUo VVOrK 


z40oooo.oUDoy 
Z40OOO / . / 4 / Do 
z404zZZ. / i Dio 
24o4zZD.71o4o 


0.00018 
0.00015 
0.00036 
0.00037 


rioiman et ai. (zuud) 
rioiman et ai. \Z\jvo) 
L-accres et al. (zUUyj 
Caccres et al. (zUUyJ 


2453808.9170 





0010 


AXA'^: Poddanv et al. f2010) 


2453832.5683 





0020 




2453879.8643 





0011 




2453887.7477 





0006 


" 


2453903.5141 





0008 


" 


2454214.8927 





0010 


" 


2454285.8403 





0010 


" 


2454553.8624 





0010 




2454620.8655 





0015 




2454620.8678 





0008 


" 


2454624.8100 





0014 




2454628.7515 





0004 




2454648.4586 





0011 




2454888.8900 





0007 




2454908.5962 





0010 


)) 


2454959.8375 





0005 




2454963.7789 





0005 


)) 


2454963.7790 





0010 


)) 


2454967.7192 





0007 




2454971.6613 





0006 




2454983.4866 





0008 




2454987.4285 





0008 




2455058.3769 





0010 




2455290.9235 





0006 




2554171.5333 





0017 


Raetz et al. f2009") 



NOTE: Times are on the BJD Terrestrial Time (TT) 
system. For times on the BJD UTC system subtract 
65.184 s for times before 2454832.5 Juhan Date (UTC) 
and 66.184 s for times after. 
^ Amateur Exoplanet Archive 

[http://brucegary.net/AXA/x.htm| 



